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ABSTRACT 



We phenomenologically put local constraints on the rotation of distant masses 
by using the planets of the solar system. First, we analytically compute the or- 
bital secular precessions induced on the motion of a test particle about a massive 
primary by a Coriolis-like force, treated as a small perturbation of first order in 
the rotation, in the case of a constant angular velocity vector directed along a 
generic direction in space. The semimajor axis a and the eccentricity e of the test 
particle do not secularly change, contrary to the inclination /, the longitude of 
the ascending node fi, the longitude of the pericenter w and the mean anomaly 
M.. Then, we compare our prediction for (vj) with the corrections Avj to the 
usual perihelion precessions of the inner planets recently estimated by fitting long 
data sets with different versions of the EPM ephemerides. We obtain as prelim- 
inary upper bounds \^ z\ < 0.0006 — 0.013 arcsec cty~^, \^x\ < 0.1 — 2.7 arcsec 
cty~^, \^y\ < 0.3 — 2.3 arcsec cty~^. Interpreted in terms of models of space-time 
involving cosmic rotation, our results are able to yield constraints on cosmolog- 
ical parameters like the cosmological constant A and the Hubble parameter Hq 
not too far from their values determined with cosmological observations and, in 
some cases, several orders of magnitude better than the constraints usually ob- 
tained so far from space-time models not involving rotation. In the case of the 
rotation of the solar system throughout the Galaxy, occurring clockwise about 
the North Galactic Pole, our results for ^ ^ are in disagreement with the expected 
value of it at more than 3 — cr level. Modeling the Oort cloud as an Einstein- 
Thirring slowly rotating massive shell inducing Coriolis-type forces inside yields 
unphysical results for its putative rotation. 

Subject headings: gravitation — celestial mechanics — astrometry — ephemerides — 
planets and satelhtes: individual (Mercury, Venus, Earth, Mars) 
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Introduction 



According to iMachI (jl883[ ). tlie fictitious forces, proportional to tlie inertial mass m 
of a body, wliich arise in a local non-inertial frame are merely due to a state of relative 
rotation ^ among the considered frame and the distant masses of the universe, contrary to 
the Newtonian picture attributing them to the rotation of the frame with respect to the 
absolute space. If at a given point the average rotation of the rest of the universe yields a 
given amount of non-inertiality, in another place the situation will be, in general, different 
because of the local average rotation of the remote masse s over there. For a c count s of 
the history of the Machian ideas , the reader is referred to iBarbour fc Pfisterl ( 1l995l ) and 
Lichtenegger &: Mashhoon ( 20071 ). Anyway, Mach neither had at disposal nor developed a 
mathematical theory supporting his views. 

The situation change d with the adv ent of the Einsteinian investigations for a relativistic 
theory of gravitation. As lM0llerl (jl969[ ) writes, Einstein advocated a new interpretation 
of the fictitious forces in accelerated systems of reference: instead of regarding them as 
an expression of a difference in principle between the fundamental equations in uniformly 
moving and accelerated systems, he considered both kinds of reference to be completely 
equivalent as regards the form of the fundamental equations: and the "fictitious" forces 
were treated as real forces on the same footing as any other forces of nature. The reason for 
the occurrence in accelerated systems of reference of such peculiar forces should, according 
to this new idea, be sought in the circumstance that the distant masses of the fixed stars are 
accelerated relative to these systems of reference. The "fictitious" forces are, thus, treated 
as a kind of gravitational force, the acceleration of the distant masses causing a "field of 
gravitation" in the reference frame considered. A recent discussion of t he rela t ivity of 
rotational motion within the context of general relativity can be found in lGr0nl (j2009[ ). For 
a detailed historical review of the consequences of rotating m asses in relati vistic physics, 
including also the centrifugal-type effects of order see iPfisterl ( 120071 ): here we will 



briefly review just some salient points. Limiting t o the larger first-order effects in rotation 
with the so-callecQ "Entwurf" tensorial theory ( Einstein fc Grossmann 1913 ). developed 
with M. Grossmann in 1913, Einstein and M. Besso considered a spherical rotating mass 
shell in the so-called Einstein-Besso manuscriptEl, pp. 36-37, and worked out the resulting 
Coriolis- l ike dr agging force occurring inside. Such a result was publicly presented by 
Einstein ( 1913 ) in his talk in September 1913 at the Naturforscherversammlung in Viennalf]. 
Einstein, presumably realizing that the Coriolis-type effects in the final version of his 
theory of gravitation would not have differed qualitatively from his previous results in the 



^It is a theory in which the left-hand side of the field equations is a "tensor" which is 
covariant only with respect to a reduced class of coordinate transformations. 

^Such an important document has been reprinted and commented in lKlein et al.l ( 119951 ). 

^H. Thirring attended it. 



-4- 



Entwurf theory (jPfisterll2007l ). did not repeat the calculations after he published his general 
relativity in November 1915. Such a task was undertaken by H. Thirring in 1917 who, 
anyway, initially neglected the ef fects of order Cp j"^) concentrating on the the much smaller 
centrifugal on es of order Oj"^^) flThirrind 119171 ). It was after the letter by Einstein of 2 
August 1917 fISchulmann et al.lll998h to him that Thirring was put on th e track o f the 
Coriolis-type dragging effects in general relativity. computed them in Thirring (1918) 
by obtaining a result only differing by a multiplicative factor of two from the Entwurf-based 
calculation by Einstein: the general relativistic dragging inside a massive rotating hollow 
shell is 

A = -2{Exv), (1) 



with 



H = g*, q 



4Gm 



3c^R' 

where G is the Newtonian constant of gravitation, dJl and R are the mass and the radius 
of the shell, respectively, and c is the speed of light in vacuum; S has the dimensions of 
s~^. The validity of eq. ([I]) carries over to all regions in which R <^ r. The adimensional 
coefficient q in eq. ([2]) was modified as 



(2) 



4e(2 



[l + e){3-e) 



with 



GM 

2^' 



(3) 



(4) 



by lBrill fc Cohen! ( 1l966l ). It is interesting to note that by choosing a suitably rotating frame 
with angular velocity it is po s sible to eli minate the Cor iolis-force inside the shell. In 
principle, the results by lEinsteiru ( 1l913l ) and iThirringI ( 1l918l ) concerning the rotating hollow 
shell can be criticized from a Machian point of view because of the asymptotic flatness of 
the exterior solution, instead of using cosmological boundary conditions. Anyway, in recent 
times it was demonstrated that, essentially, the Coriolis-like dragging effects inside the mass 
shell in an asymptotically flat back ground carry over with minor changes to cosmological 
boundary conditions. For example, iKleinl ( 119931 ) showed that it is possible to embed a 
slowly rotating massive shell with flat interior in a rotati onally per t urbed Friedmann 
universe obtaining dragging effects comparable to those by iThirringI ( 1l918l ): of course, they 
depend on the type of the Friedmann universes {k = 0, ±1), and on its mass density. 

Putting aside the localized Einstein-Thirring rotating mass shell, in the early years 
of general relativity attempts were made to construct exact or approximate solutions of 
the field equations with rotating matter source having cosmological significance; after all. 



'^T hirring did not acknowledge the contributions by Einstein to this specific topic (jPfister 
2007h . 
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rotation is an ubiquitous phenomenon in nature. We can observe rotating objects at all 
scales of the Universe, from the elementary particles to planets, stars and galaxies. The 
question is , whether this pr operty is an attribute of the whole universe at a very large scale 
structure ( IWhittakerl Il945l ) . On the other hand, if our universe does not rotate, then it 
should be explained why and how this happens. Since rotation is generic in the universe, 
the possible rotation of the universe cannot be excluded at the very beginning. Moreover, 
we should explain the physical mechan ism which prevents universal rotation of the universe 
( lObukhov et al.ll2002l ). iLanczosI (jl924j ) was the first to consider the possibility of a globally 
rotating universe modeled as a rigidly rotating dust cylinder of infinite radiu^. For a 
gene ral overview on physical foundations and observational effects of cosmic rotation, see, 
e.g., lObukhovl (|2000(). After earl y speculations concerning the role of cosmic rotation on the 



galaxy formation ( lGamowlll946l ) and the universe's structure fjvon Weizsackerlll948l ). iGodel 
( 119491 ). proposed a stationary cosmological model in which th e angu l ar vel ocity of the 



Godell ( 119491 ) considered 



cosmic rotation, is a positive constant having dimensions of s" 
the simplest matter source, i.e. ideal dust with the energy-momentum tensor T^^, 
and by solving the Einstein field equations with a cosmological constant A obtainecij 



^2 

IF 



-A. 



(5) 



Two problems of the solution by iGodell (119491 ) were the lack of expansion of the resulting 
model of the universe, and the violation of causality due to the existence of closed timelike 
curves. It was later demonstrated that suitable extensions of the original model, including, 
e.g., more gener al energy-momentum tensors or cosmic shear, may circumvent such issues 
(lQbukhovll20"00h . The cosmic rotation affects the polarization of the radiation propagating 
in the curved spacetime giving rise to an observable effect. Concerning the magnitude of 



the c osmic vorticity, it turns out to be of the order of the Hubble paramet er (ILarson et al. 



20101) Hq 
yieldB ^! 



71.0 km s-i Mpc" 
: (1.8±0.8)g n, / = 



^ = 2.3 X 10 

295° ±25°, 6 



18 



24° ± 20° on the basis o f the d ata by iBirch 



i.o mas cty-i. Indeed. lObukhovl (Il992h 



(119821 ): lObukhovl (120001 ). interpreting the data by iNodland fc RalstonI (119971 ) concerning 



the dipole effect of the rotation of the plane of polarization as arising from cosmic rotation, 
finds ^ = (6.5 ± 0.5)i/o, / = 50° ± 20°, h = -30° ± 25°. Within the errors, the directions 



of ^ are orthogonal to each other. Loca l, astronomical conseque nces of the iGodell (119491 ) 
model have been recently investigated in IWilson fc Blomd (120091 ) by looking for a possible 
explanation of the Pioneer anomaly occurring in the remote regions of the solar system. 

Recently, investigations of rotational perturbations of pure Friedmann-Lemaitre- 
Robert son- Walker (FLRW) cosmologies, by discarding also in this case the local rotating 



^A serious drawback of such a model was that the dust density diverges at radial infinity. 
^For Godel A is negative. 

^Here / and b denote the Galactic longitude and latitude, respectively. 
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2001, 2006 



mass shell, and o f their Machian impact have been produced (IBicak et al.l l2004t ISchmid 
20091). In t h e framework of the FLRW mo d els w i th gl obal rotation, 



Godlowski fc Szvdlowskil ( 120031 ): ISzvdlowski fc Godlowskil (j2003l . 120051) tackled the problem 



of the recently obser ved cosmic acceleration from sup ernovae SNIa and Cosmic Microwave 



Background (CMB). ISzydlowski fc Godlowskil (120051 ) showed that the acceleration of the 
universe can be explained in terms of the global rotation of the universe. Moreover, 
global rotation gives a natural explanation of the em pirical relation between angula r 



momentum for clusters and superclusters of galaxies (ISzydlowski &: Godlowskil l2005l ). 



Pontzen & Challinor 


2007) 


(1973 


); 


Barrow et al. ( 


1985) 



2007 ): for earlier investigations on this topic, see ICoUins fc Hawking 
19851 ) . Other cosmological features like the recent discoveries of some 
non-Gaussian properties of the Cosmic Microwave Background Anisotropies (CMB A), such 
as the suppression of the quadrupole and the alignment of some multipoles have attracted 
further attention to rotationally perturbed FLRW models with a cosmological constant 
( ISu fc Chull2009[ ) . By comp a ring t he second-order Sachs- Wolfe effect due to rotation with 
the CMBA data, ISu fc Chul (120091 ) constrained the angular speed of the rotation to be less 
than 20 mas cty~^ = 3.2 x 10~^^ s~^ at the last scattering surface. 



In all the cases considered so far, a rotation is induced with respect to a locally 
quasi-inertial Fermi frame; thus, a Corioli s-tvpe acceleratio jl affecting the motions of test 



particles moving with respect to it arises (IWilson fc Blomell2009l ). See, e.g., ISilkl ( 119661 ) for 
early studies of the Coriolis and centrifugal accelerations acting on a test particle moving 
in the Godel spacetime. This fact opens, in principle, interesting perspectives to put local 
constraints on the angular velocity vector in a purely phenomenological way from local 
dynamics of the planets of the solar system. Recall that the celestial reference system (ICRS) 
used is based on a kinematical definition, making the axis directions fixed with respect to 
the distant matter of the universe. The system is materialized by a celestial reference frame 
(ICRF) defined by the precise coordinates of extragalctic objects, mostly quasars, BL Lac 
sources and few active galactic nuclei (AGNs). The current positions are known to better 
than a milliarcsecond (mas)= 4.8 x 10~^ rad. According to the lAU recommendations, the 
origin is to be at the barycenter of the solar system, and the directions of the axes should 
be fixed with respect to the quasars via VLBI observations. Such recommendations further 
stipulate that the principal plane should be as close as possible to the mean equator at 
J2000.0 and that the origin of it should be as close as possible to the dynamical equinox of 
J2000.0. It turns out that the uncertainty from the representation of the ICRS is smaller 
than 0.01 mas, and the axes are stable to ±0.02 mas. Note that this frame stability is 
based upon the assumption that the extragalctic sources have no proper motion and that 



^This is a first-order effect, not to be confused with the second-order ones due to the 
tidal forces induced in th e local frame by the background cosmological space-time; see, e.g.. 



Cooperstock et al.l (119981 ). 
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there is no global rotation of the universe f lMc Carthy fc Petitll2004l ). It must be noted that 
the aforementioned accuracies make, in principle, meaningful the analysis proposed; indeed. 



^ = V-Ac = 3.4 X 10 
cosmological constant 



-18 



2.2 mas cty i n view of the current ly accepted value of the 



-A = 1.26 X 10-^2 m-2 (iHinshaw et al.ll2009l ). 



The paper is organized as follows. In Section [2] we analytically work out the 
long-term effects on the orbital motion of a test particle perturbed by a small Coriolis-type 
acceleration. In Section [3] we compare the resulting predictions with the latest observational 
determinations of the non-standard perihelion precessions of some planets of the solar 
system. Section H] is devoted to summarizing our results and to the conclusions. 



2. Analytical calculation 



The orbital effects induced on a test particle by any acceleration A quite smaller 
than the dominant Newtonian monopole —GM/r^ of the primary of mass M, whatever its 
physical origin may be, can be worked out with standard perturbative techni ques by using, 
e.g., the Gauss equations for the variations of the Keplerian orbital elements flBertotti et al. 



20031). They are 



da 
H 



2 

nrj 



eAnsmf + AT ( 



p 
r / J 



^ = ^{Ansinf + AT 
dt na 



(6) 



(7) 



dl 

lu 

dn 
Itt 

duj 

m 

dM 

IT 



An (- ] cosM, 

nar] 



Am ( - 1 sinw 



na sin It] 
V 



nae 



-Ar cos / + At ( 1 + - ) sin / 



— cos/ 



dn 
11' 



n Ar - 

na \a 



/ doj ^dQ\ 



(9) 
(10) 

(11) 



where a,e,I,Q,uj,Ai are the semimajor axis, the eccentricity, the inclination, the longitude 
of the ascending node, the argument of pericenter and the mean anomaly, respectively, of 
the orbit of the test particle. Moreover, / is its true anomaly reckoned from the pericentre 
position, M = a; -|- / is the argument of latitude, n = \/ GMj a?" = 27r/Pb is the unperturbed 
Keplerian mean motion related to the unperturbed Keplerian orbital period P^, r) = a/1 — 
and p = a(l — e^) is the semi-latus rectum. Finally, Ar, A^, A^ are the projections of the 
perturbing acceleration A onto the radial R, transverse T and out-of-plane N directions of 
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the particle's co-moving frame {f, f, u}. Since in the following we will work out the net, 
secular effects of A on the particle's Keplerian orbital elements, the right-hand sides of eq. 
(EJ-eq. ( ITT]) have to be evaluated onto the unperturbed Keplerian trajectory. To this aim, it 
will turn out to be convenient to use the eccentric anomaly E instead of the true anomaly 
/; useful conversion relations are 

r = a(l-ecosE), (12) 
cosE-e 

cos/ = (13) 

1 — e cos E 

sm/ = (14) 

1 — e cos E 

(1 — ecosE) ,^ . ^, 

dt = ^ ^-dE. 15 

n 

In our case, the disturbing acceleration is a phenomenological Coriolis-like one 

A = -2* X V, (16) 

where \Ef has the dimensions of s~^, and, in general, we do not make any a-priori assumption 
on the rotation velocity vector i.e., we pose 

* = + ^yj + <if,k. (17) 

Actually, in some models of rotating universe depends both on time and radius (jSu fc 



20091 ). but given the typical temporal and spatial scales of the solar system orbital motions, 
we can safely assume it to be a constant vector. The components of the planet's velocity v 
entering eq. ( fT6l) can be evaluated onto the unperturbed Keplerian ellipse as 

dx dE . . 

= dEW ^''^ 

dz dE . . 



in which 



X = r(cosf2cosn — cos / sin sin u) , (21) 
y = r (sin f2 cos M + cos / cos sin m) , (22) 
z = r sin/ sin M, (23) 
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and 



dE 



n 



(24) 

1 - e cos E ^ ' 

Note that, since we are going to take an average over one orbital revolution of the right-hand 
sides of the Gauss equations, we can consider 7,^2, w as constant, so that the partial 
derivatives of x, z with respect to E involve only r and u. 

The time-varying unit vectors of the co-movin g frame along the radial, transverse and 
normal directions are ( iMontenbruck and Gillll2000l ) 



cos VL cos u — cos / sin VL sin u 
sin VL cos u + cos I cos VL sin u 
sin I sin u 



(25) 



sin u cos Vt — cos / sin Vt cos u 
sin Vt sin u + cos / cos Vt cos u 
sin J cos u 



V 



sin / sin f2 
- sin / cos 
cos / 



Thus, the projections of the perturbing acceleration on them, defined as 





= A- 


f 


At 


= A- 


f 


An 


= A- 





are of the form 



Ai= Yl Vijia,e,I,n,u;E)^j, t = R,T,N, 

j={x,y,z} 

where the coefficients Vij have dimensions of velocities. The coefficients of the radial 
acceleration An are 



(26) 



(27) 



(28) 
(29) 
(30) 

(31) 



V 



Rx 



V 



Ry 



V 



Rz 



2o?nri sin / sin VL 



2o?nri sin / cos VL 



2a?nri cos / 



(32) 
(33) 
(34) 
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The transverse acceleration At is built of 



2a'^ensinl sinQsinE . 

vtx = , (35) 

r 

^, 2a^e?7, sin / cos r2 sin 

Vry = , 36) 

r 

2a^er2COS Jsini? . , 

Vtz = . (37) 

r 



The coefficients of the normal acceleration are 

^Nx = [sin£' (cosfisincj + cos/sinficoso;) — 

r 

—r] cos E (cos Q cos u — cos / sin Q sin u)] , (38) 

Vat^ = [sin E (cos / cos Q cos a; — sin Q sin cj) + 

r 

+?7 cos -E (cos / cos sin a; + sin Q cos u)] (39) 
^Qj^Ti sin 7^ 

Vatz = (sinii^coso; + r/cosE'sino;) . (40) 

In order to make contact with the latest observational determinations from planetary 
motions, it is convenient to work out the secular precession of the longitude of the pericenter 
zu defined as 

zu = n + u. (41) 

From eq. dHJ-eq. f IlOp it turns out that, actually, a small perturbing Coriolis-like acceleration 
induces a non-zero secular precession of w given by 

(^)= E ^.■(^'^)*.' (42) 

j=x,y,z 

with the adimensional coefficients Vj given by 

= - tan sin Q, (43) 

Vy = tan cos fi, (44) 
= -1. (45) 
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Also other Keplerian orbital elements, for which no observational investigations yet 

exist, undergo secular changes. It turns out that no secular variations occur for the 
semimajor axis a and the eccentricity e. Instead, the inclination I experiences a secular 
rate given by 

(/) = ^.(^)^.' (46) 

with 

Zj. = — cos (47) 

Xy = -sinn, (48) 

X. = 0. (49) 

E ^^•(^'^)^^' (50) 

Mx = cot/ sin f2, (51) 

My = -cot/cosfi, (52) 

K = -1. (53) 
The secular change of the mean anomaly Ai is 

'M)=n+ J2 -^j(^'^)*j' (54) 



The node fl precesses at a rate 



with 



with 



Aix = — ?7 CSC / (cos / + 3 sin^ /) sinf2, (55) 
M.y = ?7 (cot / + 3 sin /) cos (56) 
M, = r/ (1-3 cos J). (57) 
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Confrontation with the latest observational determinations 



Recently, iPitjeval f l2010l ) has analyzed more than 550000 planetary observations of 
several kinds covering the time interval 1913 — 2008. She used the dynamical force models 
of the EPM2008 ephemerides by estimating about 260 parameters. The reference frame 
used to numerically integrate the equations of motion, assumed non-rotating with respect 
to extra-galactic sources and, thus, locally inertial, is a barycentric one aligned with ICRF 
by including into the total solution also the VLBI data of spacecraft near the planets; 
actually, the fundamental plane is the mean ecliptic at J2000.0, so that the z axis points 
towards the North Ecliptic Pole, in Draco constellation, with right ascension a = 18*^0^0. 0*^ 
and declination 6 = +66°33 38.6 . Among the estimated parameters, she also determined 
corrections Ati; to the standard Newtonian/Einsteinain perihelion precessions of all the 
planets of the solar system including Pluto as well; such corrections, by construction, 
account for any unmodelled/mismodelled dynamical effects, so that they can be used, in 
principle, to preliminarily put constraints on ^. Actually, the entire planetary data set 
should be re-processed with ad — hoc modified models to account for the Coriolis-type 
effect investigated; one or more dedicated solve-for parameters should be simultaneously 
estimated along with all the other ones, but this lies outside the scopes of the present 
paper. In order to constrain we will use the inner planets, whose estimated perihelion 
corrections are listed in Table [H because they are more accurate. In Table [3] we quote 
the coefficients Vj of the Coriolis-type perihelion precessions for the inner planets whose 
relevant Keplerian orbital elements are listed in Table [2J 

By equating the corrections Atu of Table [T] to the predicted Coriolis-type precessions 
(w) of eq. (H2|) for three planets it is possible to set up a non-homogeneous algebraic linear 
system in the three unknowns '^xj'^y,'^ z- By using Mercury, Venus and the Earth allows 
to obtain 

= -0.9 ± 1.8 arcsec cty-\ (58) 
= -1.3 ±2.0 arcsec cty"\ (59) 
= -0.006 ± 0.007 arcsec cty-\ (60) 



Table 1: Estimated corrections Aw, in mas cty~^, to the standard perihelion precessions 
with the EPM2008 ephemeri des. The quote d errors are not the formal, statistical ones but 
are realistic. From Table 8 of iPitjeva fl2010h . 



Mercury 


Venus 


Earth 


Mars 


-4 ±5 


24 ±33 


6±7 


-7 ±7 
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with 



\1/ = 1.6 ± 2.6 arcsec cty 



(61) 



The quoted uncertainties have been obtained by hnearly propagating the errors in Aw of 
Table [T] It turns out that the other combinations involving Mercury, Venus and Mars, and 
Venus, Earth and Mars yield similar constraints, while those from Mercury, Earth and Mars 
are about one order of magnitude weaker. The tightest constrain holds for {"^^l < 0.013 
arcsec cty"^; instead, may be as large as ~ 1 arcsec cty~^. Recall that 1 arcsec cty~^ 



1.5 X 10 



-15 



In Table H] we quote the values of the corrections Aw estimated by iPitjeval ( l2005l ) with 
older versions of the EPM ephemerides and less extended data sets; for Venus, Earth and 
Mars they are, for some reasons, more accurate by about one order of magnitude than the 
more recent results of Table [1] The data of Table HI applied to Mercury, Venus, Earth, yield 



-1- ri 



-0.02 ± 0.08 arcsec cty^ 



(62) 



-0.1 ± 0.2 arcsec cty" 



(63) 



corresponding to 



= -0.0002 ± 0.0004 arcsec cty" 



= 0.1 ± 0.2 arcsec cty \ 



(64) 



(65) 



of the order of 10 s ^, i.e. one order of magnitude more accurate that the limit from eq. 

(EH). 



By comparison, iFolknerl ( 120101 ) phenomenologically tested if the locally inertial 
reference frame in which the orbits of the planets are usually integrated is not rotating 
with respect to the rest of the universe. He compared the mean motion of Mars relative 
to Earth determined from ranging measurements with the mean motion determined from 
VLBI o bservations of Mars-orbiting spacecraft relative to stars or extra-galactic radio 
sources. iFolknerl ( l2010l ) finds that the dynamical rotation rate of the solar system relative to 



Table 2: Longitude of the ascending node f2 and inclination /, both in deg, for the inner 
planets. Reference frame: ICRF/J20 00. Coordinate systein : Ecli ptic and Mean Equinox of 
Reference Epoch. From Table A. 2 of iMurray and DermottI (119991 ). 





n 


/ 


Mercury 


48.33167 


7.00487 


Venus 


76.68069 


3.39471 


Earth 


348.73936 


0.00005 


Mars 


49.57854 


1.85061 
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extra-galactic radio sources is less t han 0.00 4 arcse c cty~^ = 6 x 10~^^ s~^. Actually, more 
details on the approach followed by iFolknerl ( l2010l ) would be needed to make a meaningful 
comparison with our results. In fact, no long-term effects may occur for the mean motion 
because the Coriolis force does not secularly affect the semimajor axis a; instead, the mean 
longitu de £ = zu + Ai experiences a secular variation. A result quantitatively similar to 
that by iFolknerl (120101 ) could be obtained with our method by a-priori setting = \E'y = 



because, in this case, it would simply be 



(66) 



-3cos/^. 



(67) 



Interpreted in terms of rotation of the universe (jWilson &: Blomd 120091 ). our results 
show that local astronomical effects in the solar system do not show statistically significant 
evidence for its existence, putting independent upper bounds on its magnitude. They are 
not in contrast with, e.g., the Godel- type expec t ed va lu e \1' = 2 rnas cty ~^ an d with th e 
upper bound of 0.02 arcsec cty"^ by ISu fc Chul fcoOOh . lObukhovl (Il992[ ) and lObukhov 



( I2OOOI ) yield not only the magnitude of the cosmic rotation but also its spatial direction. 
More precisely, from the Galactic longitude / and latitude fe of ^' released by lObukhov 
( 119921 ) it is possible to infer its direction with respect to the solar system barycentric frame 
used 



-0.73 ±0.24, 



(68) 



-0.51 ±0.34, 



-0.45 ±0.32, 



(69) 
(70) 



Table 3: Coefficients Vj, j = x,y,z for the inner planets computed from Table [2] and eq. 
gSD-eq. (USD. 





V. 




V, 


Mercury 


-0.04572069 


0.04069034 


-1 


Venus 


-0.02883601 


0.00682681 


-1 


Earth 


8 X 10-s 


4.3 X 10-^ 


-1 


Mars 


-0.01229570 


0.01047240 


-1 
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so that, by summing in quadrature the errors in \1/ and in the components of the unit vector, 

= -1.97 ± 1.09 mas cty"\ (71) 

= -1.38 ± 1.10 mas cty"\ (72) 

= -1.21 ± 1.02 mas cty'^ (73) 
It can be noted that eq. f l7T]) -eq. f l75]) are compatible wit h both eg. fl^-e a. fIBCT]) and the 



more stringent bounds of eq. (!62|) -eq. ( 164|) . The values by lObukhovl ( l2000l ) for / and b yield 
for the components of the unit vector of ^ 

= 0.72 ±0.27, (74) 
= -0.63 ±0.30, (75) 



Thus 



= 0.27 ±0.34. (76) 

= 7.0 ± 2.7 mas cty"\ (77) 
= -6.1 ± 2.9 mas cty-\ (78) 

= 2.6 ± 3.3 mas cty~\ (79) 

Also in this case, eq. fl77|) -eq. f l79|) are compatible with both eq. f l58|) -eq. f l60|) and 
eq. fl5^ -eq. f lM|l . It is remarkable to note that our approach, in the framework of the 
polarization radiation due to cosmic rotation, is able to put local, solar system-scale 
constraints on the Hubble parameter Hq which are just 2 — 3 orders of magnitude larger 
that its value determined from cosmological observations. Incidentally, let us note that 
the relation between the Godelian universe's vorticity and the cosmological constant 



Table 4: Corrections Aw, in milliarcsec cty ^, to the standard perihelion precessions esti- 
mated by E.V. Pitjeva with the EPM2005 (Mercury, Earth, Mars) and EPM2006 (Venus) 
ephemeri des. The quote d errors are not the formal, statistical ones but are realis tic. From 
Table 3 of Eltjeval (|2005|) (Mercury, Earth, Mars) and Table 4 of IPienga et"all fl2010l ) (Venus). 



Mercury 


Venus 


Earth 


Mars 


-3.6 ±5.0 


-0.4 ±0.5 


-0.2 ±0.4 


0.1 ±0.5 
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allows us to put local constraints on it much tighter than those usually obtained from the 
Schwarzschild-de Sitter metric usually adopted so far. Indeed, from eq. ( !6T]) it follows 



A < 2 X 10 



-46 -2 

m , 



while eq. ( l65|) yields 



- A < 1 X 10 



-48 



m 



(80) 



^1) 



on the contrary, the solar system-based constraints obtained from the perturbations by the 
Hoo ke-type radial acceleration due to the Schwarzschild-de Sitter metric are of the or der 
o^ flloriol l: 



2006 : iJetzer fc Serendl2006t iKagramanova et al.ll2006t ISereno fc Jetzerl 120061 



A < 10 



-40 



10-^2 m- 



^2) 



With regard to the connection with cosmological parameters, interesting perspectives 
may open up in near-mid future if and when the interplanetary laser ranging technique 
( iDegnanI 120081 ) will be implemented allowing for a notable improvement of the planets'orbit 
determination. 

Concerning the application of our Coriolis-type res ults to the rotation of the solar 
system through the Milky Way (jWilson fc Blomd 120091). let us start by n oting that it 
occurs clockwise about the North Ecliptic Pole, (IBinney fc Tremaind 119941 ). located in the 
constellation Coma Berenices. Thus, the associated unit vectoi0 4f has ecliptic longitude 
and latitude 



A 



0.02319 deg. 



(83) 



and components 



= -29.81149 deg 

= 0.867, 
= 3 X 10-^ 
= -0.497. 



(84) 

(85) 
(86) 
(87) 



^Here A is defined positive. 

^^Conventionally, it is directed towards the South Galactic Pole, located in the constella- 
tion Sculptor, so that the Sun's rotation appears anticlockwise from its tip. 
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The magnitude of the Galactic angular velocity can approximatelj0 be evaluated as 

Vq/ {27rrQ), so that 



[IM ± 0.19) X 10"^^ s"^ = 0.101 ± 0.009 arcsec cty"\ 



(88) 



where we used Vq = (254 ± 16) km s'^, Vq = (8.4 ± 0.6) kpc (IReid et al.ll2009l ). Its 
components are 



0.087 ±0.008 arcsec ctj-\ 



(89) 



C(10-^) arcsec cty-\ 



(90) 



-0.050 ± 0.004 arcsec cty-\ 



(91) 



It must be noted that eq. fl9T]) is statistically incompatible with both eq. (160!) and eq. (E 
at more than 3 — 0" level. 

It may be interesting to point out that if we interpret our results in terms of an 



Einstein-Thirring rotating massive shell identified wit h the Port cloud (jOpiklll932l : lOort 



l-Oort 



38mffi and 



1950l ). we get unphysical results. Indeed, by assuming flWeissmanlll996l ) TIq 
R ^ 10^ au, eq. yields 

goort ~ 1.5 X 10-1^ (92) 

Thus, eq. and eq. (^5^ would yield an angular velocity of the Oort shell as large 
as 1 — 10 s~^. Of course, it is very daring to identify an extended object of certainly 
non-uniform density like the Oort cloud with the Thirring's infinitely thin massive shell of 
uniform density. 



4. Summary and conclusions 

We analytically worked out the secular precessions of the Keplerian orbital elements of 
a test particle affected by a small extra-acceleration of Coriolis type treated perturbatively 
with the Gauss variational equations. It could be due to a relative rotation of the reference 
frame considered with respect to distant masses; we reviewed several theoretical local and 
cosmological scenarios leading to such an effect. In the calculation we did not make any 
a-priori assumptions on the spatial orientation of the angular velocity vector ^ which was, 
thus, treated as a constant vector directed along a generic direction in space. We found 
that the semimajor axis a and the eccentricity e do not undergo secular changes, while 
the inclination J, the longitude of the ascending node Q, the longitude of the pericenter cjj 



A circular motion is assumed. 
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and the mean anomaly Ai experience secular variations which are linear combinations the 
components of ^ with coefficients depending on e, /, Q. 

Then, we compared the Coriolis-induced theoretical prediction of the precession of 
w with the latest observational determinations of the corrections Aw to the standard 
Newtonian-Einsteinian precessions of the inner planets of the solar system estimated by 
fitting long data sets with different versions of the EPM ephemerides. No statistically 
significative evidence for a non-zero angular velocity ^' was found; thus, we were able 
to constrain ^. The tightest bounds occur for whose magnitude can be as large as 
0.0006 — 0.013 arcsec cty~^. The other components can be as large as {"^^l < 0.1 — 2.7 
arcsec cty~^ and < 0.3 — 2.3 arcsec cty~^. It must be pointed out that such upper 
bounds should be considered just as order-of-magnitude figures; indeed, the entire planetary 
data set should be re-processed by explicitly modeling the effect we are interested in, and 
dedicated solve-for parameters should be simultaneously estimated in a new global solution 
along with all the other ones routinely determined. 

It turns out that several values of the cosmic rotation obtained from cosmological 
observations interpreted in terms of different theoretical models are compatible with our 
results. In the framework of the rotating Godel universe, they yield local, astronomical 
constraints on the cosmological constant A of the order of 10~^^ — 10~^* m~^, i.e. several 
orders of magnitude tighter than those usually obtained so far from the Schwarzschild-de 
Sitter metric. On the other hand, interpreting our results in terms of the polarization of 
radiation propagating in a rotating universe allows to put solar system-scale constraints 
on the Hubble parameter Hq which are just 2 — 3 orders of magnitude larger that its 
value derived from cosmological observations. This opens interesting perspectives in 
view of the expected future improvements in the planetary orbit determination from the 
implementation of the interplanetary laser ranging technique. Concerning the effect of the 
rotation of the solar system through the Milky Way, occurring clockwise about the North 
Galactic Pole, it turns out that the corresponding component along the z axis of \l/ is in 
disagreement with our values for at more than 3 — a level. Finally, let us mention that 
the application of our investigations to the Oort cloud, modeled as an Einstein-Thirring 
rotating massive shell inducing Coriolis-type forces inside, yields unphysical values for its 
putative rotation. 
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